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PREFACE 


The  accurate  prediction  of  unsteady  air  loads  is  essential  to  avoiding  problems  and 
assuring  safety  in  the  many  interdisciplinary  regions  involving  aeroelasticity  and  the 
dynamics  of  active  controls.  However,  the  methods  to  predict  such  airloads  are  complex 
and  intricate.  It  was  considered  essential  to  specify  standard  configurations  and  para- 
meters, and  to  encourage  pioneering  NATO  scientists  to  report  their  results  early.  This 
would  provide  bases  for  evaluation  and  improvement  of  later  and  following  developments 
by  other  countries.  The  first  cooperative  effort  involved  isolated  surfaces  in  subsonic  ana 
supersonic  flow,  and  is  reported  by  D.L.Woodcock  in  AGARD  Report  No.583  (“A 
Comparison  of  Methods  Used  in  Lifting  Surface  Theory”,  1971).  The  noticeable  success 
led  to  another  effort  on  interfering  lifting  surfaces  and  is  reported  by  W.P.Rodden  in 
AGARD  Report  No.643  (“A  Comparison  of  Methods  Used  in  Interfering  Lifting  Surface 
Theory”,  1976).  Both  are  supplements  to  the  AGARD  Manual  on  Aeroelasticity,  Vol.Vl. 
The  latter  effort  and  report  has  also  proved  to  be  highly  successful. 


New  developments  are  rapidly  emerging  in  unsteady  aerodynamics.  The  aeroelastician 
will  continue  heavy  reliance  on  prediction  of  airloads  from  theoretical  methods.  A terse 
description  of  the  new  state  of  the  art  was  required  and  has  been  very  competently  pro- 
vided by  Dr  Rodden  in  this  report.  The  Sub-Committee  on  Aeroelasticity  and  Unsteady 
Aerodyanmics  will  use  this  framework  to  establish  high  priority,  cooperative,  comparative 
computational  AGARD  programs.  Transonic  and  control  surface  unsteady  aerodynamics 
are  likely  selection  candidates. 

fa  p. 

W.J.MYKYTOW 
Former  Chairman 
Sub-Committee  on  Aeroelasticity 


and  Unsteady  Aerodynamics 
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STATE-OF-THE-ART  IN  UNSTEADY  AERODYNAMICS 

William  P.  Rodden 
Consulting  Engineer 
255  Starlight  Crest  Drive 
La  Canada,  California  91011 
U.S.A. 


SUMMARY 

A brief  survey  of  new  developments  in  unsteady  aero- 
dynamics is  made  as  a proposal  for  establishing  another 
comparative  computational  AGARD  program.  Candidate  topics 
include  supersonic  interference,  transonic  flow,  wing-body 
interference,  control  surfaces,  rotary  loads  on  T-tails, 
interference  effects  of  vortex  shedding,  and  rotating  blades. 

A selected  bibliography  is  presented  for  each  topic  to 
illustrate  the  present  state-of-the-art  and  its  near-term 
future  potential. 

INTRODUCTION 

This  paper  is  a survey  of  the  state-of-the  art  of  unsteady  aerodynamics  only  in  a 
limited  sense.  The  survey  is  not  intended  to  be  a critical  evaluation  of  the  references 
that  constitute  the  current  status,  but  rather  to  outline  areas  which  may  be  of  general 
interest  to  the  Structures  and  Materials  Panel  so  that  another  comparative  computational 
study  might  be  made  similar  to  those  of  Woodcock1  and  Rodden2  on  comparisons  of  unsteady 
aerodynamic  calculations  for  isolated  and  interfering  lifting  surfaces,  respectively. 

Recent  surveys  of  research  problem  areas  have  been  presented  by  Laschka3  and 
McCroskey4 . These  discuss  more  advanced  topics  than  those  mentioned  here. 
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Seven  problem  areas  are  proposed  here  and  related  recent  papers  are  cited.  The 
first  six  areas  concern  fixed  wing  aircraft  and  the  seventh  considers  the  effects  of 
rotation  in  rotorcraft  and  turbomachines.  The  areas  suggested  for  fixed  wing  aircraft 
include  further  supersonic  studies,  renewed  transonic  studies,  wing-body  interference, 
control  surfaces,  rotary  loads  on  T-tails,  and  interference  effects  of  leading-  or  side- 
edge  vortex  shedding. 


SUPERSONIC  INTERFERENCE 

This  subject  was  one  of  the  comparisons  in  Ref.  2 but  only  limited  calculations  by 
two  methods  were  compared.  A large  number  of  new  approaches  ha"e  appeared  recently  and 
it  would  seem  reasonable  to  invite  applications  of  these  new  methods  to  some  of  the 
standard  AGARD  configurations.  The  new  methods  with  which  the  author  is  familiar  include: 
the  integrated  potential  formulation  of  Appa  and  Jones5,  the  kernel  function  method  of 
Cunningham* , the  finite  element  doublet  method  of  Giesing  and  KSlmdn7,  the  finite  element 
method  of  Morino,  et  al.8'9,  the  doublet-lattice  method  of  Brock  and  Griffin10,  the  con- 
stant pressure  panel  method  of  Roger11,  the  modified  Mach  Box  method  of  Chipman12,  and  the 
characteristic  box  method  of  Schmid13.  There  are  significant  differences  among  these 
new  approaches  and  the  extent  to  which  computed  results  by  each  agreed  would  add  con- 
siderably to  an  improved  understanding  of  supersonic  interference. 

TRANSONIC  FLOW 

Sonic  calculations  by  a collocation  method  were  presented  in  Ref.  1 for  isolated 
wings.  Much  progress  in  unsteady  transonic  flow  theory  has  occured  since  the  calcula- 
tions in  Ref.  1 were  collected.  A recent  survey  by  Tijdeman14  has  summarized  the  sig- 
nificant developments  of  the  past  three  years.  Table  1 of  Ref.  14  is  reproduced  here 
as  Table  1 and  the  20  references  are  also  cited  here15-34.  Most  of  these  works  concern 
two-dimensional  airfoils  but  a few  treat  the  three-dimensional  case. 

In  his  concluding  remarks,  Tijdeman  notes  "that  inserting  the  effects  of  thickness, 
incidence  of  the  airfoil  and  of  shock  waves  in  inviscid  theory  is  an  important  step 
forward.  However,  as  shown...,  the  effect  of  the  boundary  layer  is  of  the  same  order  of 
magnitude  and  thus  of  the  same  importance.  Therefore,  realistic  improvements  from  the 
aeroelastician  point  of  view  may  be  expected  only  if  also  the  second  step  is  made, 
namely  the  inclusion  of  the  boundary  layer.  In  this  respect  the  study  of  the  behavior 
of  boundary  layers  under  unsteady  flow  conditions  will  become  important."  The  importance 
of  viscous  effects  has  also  been  observed  in  hinge  moment  predictions  for  control  sur- 
faces by  Gray  and  Davies35,  and  will  be  discussed  below,  but  it  is  interesting  to  note 
here  that  thickness,  incidence  (which  may  have  larger  effects  than  thickness),  and 
viscous  effects  should  be  treated  together  (or  all  three  effects  neglected)  in  achieving 
experimental  correlations. 

The  variety  of  new  transonic  techniques  can  also  be  evaluated  by  establishing  new 
standardized  configurations.  Two-  and  three-dimensional  configurations  might  be  based 
on  the  NACA  64A410  airfoil  analyzed  so  thoroughly  by  Magnus  and  Yoshihara2 * ' 2 * , or  the 
NLR  7301  airfoil,  which  is  representative  of  the  new  generation  of  supercritical  airfoils, 
and  is  currently  undergoing  extensive  testing  and  analysis.  Reynolds  numbers  as  well  as 
Mach  numbers  should  also  be  selected  among  the  governing  parameters. 
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WING-BODY  INTERFERENCE 

A survey  of  wing-body  aerodynamic  interaction  was  presented  by  Ashley  and  Rodden35 
in  1972.  A number  of  advances  at  subsonic  speeds  have  been  made  since  then.  An  inter- 
ference theory  has  been  developed  by  Giesing  et  al.37  in  which  the  lifting  surfaces  are 
treated  by  a combination  of  the  Doublet-Lattice  Method,  Slender  Body  Theory,  and  the 
Method  of  Images.  A panel  method  for  steady  flow  for  general  configurations  has  been 
developed  by  Labrujere  and  Sytsma3'.  Bennekers,  Roos,  and  Zwaan 5 ’ have  combined  the 
Doublet-Lattice  Method  with  an  unsteady  source  method  for  bodies  to  obtain  a general 
method  for  oscillating  wing/store  configurations.  Further  applications  of  this  new 
procedure  to  oscillatinq  wing/body  combinations  have  been  made  by  Roos,  et  al1'0.  The 
method  of  Morino  et  al.  '*  can  also  be  applied  to  wing-body  combinations  at  both  subsonic 
and  supersonic  speeds. 

A standard  configuration  for  wing-body  interference  calculations  should  include  a 
fuselage  and  an  optional  store  (nacelle) /pylon  combination. 

CONTROL  SURFACES 

The  accuracy  with  which  potential  theory  has  been  observed  to  predict  the  loading 
on  interfering  lifting  surfaces  or  wing-body  combinations  decreases  toward  the  trailing 
edge  and  the  combined  effects  of  thickness  and  viscosity  (see  Tijdeman’s  remarks  above 
under  Transonic  Flow)  probably  explain  discrepancies  with  experiment  in  hinge  moment 
predictions.  Accurate  hinge  moments  are  necessary  to  determine  power  requirements  in 
active  control  systems  for  load  alleviation  and  flutter  suppression,  as  well  as  to 
predict  flutter  involving  control  surface  motion.  Some  of  the  references  cited  in  the 
preceding  section  account  for  thickness  of  wings  as  well  as  bodies.  The  additional 
effects  of  viscosity  via  a shear  layer  have  been  considered  by  Dowell,  Ventres,  and 
Yates" 1 ' * 2 ' * 3 . The  shear  layer  is  an  inviscid  approximation  to  a viscous  boundary  layer. 

ROTARY  LOADS  ON  T-TAILS 

The  rotary  loads  are  as  important  in  predicting  T-tail  flutter  characteristics  as 
the  rotary  stability  derivatives  are  in  predicting  lateral-directional  stability  and 
control  characteristics  of  aircraft.  Some  of  these  require  the  ability  to  predict  the 
spanwise  distribution  of  induced  drag  of  oscillating  surfaces.  Others  require  the 
rolling  moments  caused  by  oscillatory  yaw  and  sideslip.  The  yawing  moment  induced  by 
steady  roll  was  discussed  briefly  by  K£lm£n  et  al.""  and  Hancock"3  . An  accurate  method 
for  predicting  spanwise  distribution  of  induced  drag  in  steady  subsonic  flow  has  been 
developed  by  Lan"*,  and  has  been  extended  to  the  case  of  oscillatory  motion  in  two-*7 
and  three-dimensions*’.  The  oscillatory  three-dimensional  case  permits  the  analysis  of 
ornithopter  propulsion*9  (as  well  as  that  of  fish  and  birds!)  as  well  as  the  oscillatory 
rotary  loads  on  a T-tail.  The  oscillatory  propulsion  of  interfering  lifting  surfaces  in 
two-dimensional  flow  has  also  been  considered  by  Bosch50. 

The  problem  of  spanwise  loading  and  rolling  moments  caused  by  oscillatory  yawing 
and  sideslipping  has  been  studied  by  Isogai  and  Ichikawa51,  but  only  for  incompressible 
flow. 


The  rotary  aspects  of  T-tails  have  not  been  studied  at  Mach  numbers  above  subsonic. 

The  standard  AGARD  T-tail  configuration  can  provide  a basis  for  comparison  of  cal- 
culated rotary  loads. 

INTERFERENCE  EFFECTS  OF  VORTEX  SHEDDING 

Recent  interest  has  been  shown  regarding  the  effects  on  steady  loads  and  flutter  of 
the  leading-edge  vortex  separation  from  a highly-swept  leading  edge  of  either  the  delta 
or  the  swing-wing  configuration.  This  is  a nonlinear  problem  in  the  flow  field  geometry. 
Recent  progress  in  predicting  steady  loads  at  subsonic  speeds  by  panel  and  lattice 
methods  has  been  made  by  Brune,  et  al.52  and  by  Kandil,  et  al.53  and  a very  complete  bib- 
liography is  contained  in  Ref.  53.  The  unsteady  subsonic  case  is  currently  under  investi- 
gation. 

The  interference  between  wing-tip  vortices  and/or  wing  jet-flap  wakes  and  horizontal 
tail  load  distributions  has  also  attracted  recent  attention.  The  steady  flow  problem 
at  subsonic  speeds  has  been  investigated  by  Goldhammer  et  al.5*  and  by  Shollenberger55 . 
Extensions  to  unsteady  flow  are  also  under  investigation55'57.  Becker's  investigation5 7 
of  slender  wing-tail  configurations  gave  good  correlation  between  prediction  and  measure- 
ment, but  the  predicted  airforces  were  usually  high. 

The  existing  standard  AGARD  configurations  are  adequate  for  comparison  of  calculated 
nonlinear  wake  interference  effects,  although  a horizontal  tail  with  a smaller  span  than 
the  wing  might  be  a more  practical  configuration. 
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ROTATING  BLADES 

Significant  progress  has  been  made  in  adapting  lifting  surface  theories  to  predict 
the  loading  of  helicopter  rotors.  References  27,  28,  58-63  present  various  approaches  and 
applications  that  all  show  promise  as  a replacement  for  Strip  Theory.  Numerical  comparisons 
would  appear  to  be  timely  in  view  of  the  importance  of  the  helicopter  in  V-STOL  Technology. 

The  Proceedings  of  a Workshop  on  Aeroelasticity  in  Turbomachines6 4 indicates  that 
the  cascade  problem  is  only  being  analyzed  by  two-dimensional  methods.  Recent  develop- 
ments in  subsonic  cascade  theory  for  staggered  compressor  rotor  blades  have  been  made 
by  Rao  and  Jones66  and  by  Jones  and  Moore66.  New  solutions  for  supersonic  cascades  have 
been  obtained  by  Verdon  and  McCune67  and  by  Yates6*. 

CONCLUDING  REMARKS 

A brief  survey  has  been  presented  that  illustrates  a number  of  new  aspects  and 
refinements  that  have  been  made  recently  in  unsteady  aerodynamics  for  interfering  con- 
figurations. Each  topic  is  not  only  an  interesting  problem  in  its  own  right,  but  also 
is  an  important  part  of  the  whole  problem  of  aerodynamic  configurations,  and  the 
analysis  of  each  has  achieved  a reasonable  level  of  sophistication.  This  level  of 
sophistication  suggests  that  a program  of  comparing  calculated  results  is  feasible  and 
could  be  a profitible  experience  for  the  contributors  and  beneficial  to  members  of  AGARD, 
in  the  same  way  that  Refs.  1 and  2 have  been  useful. 

Configurations  and  parameters  would  have  to  be  agreed  upon.  A further  study  of 
supersonic  interference  among  lifting  surfaces  requires  no  new  configurations.  Tran- 
sonic flow  studies  would  require  a definition  of  thickness  distributions;  the  NLR  7301 
airfoil  would  be  a good  choice.  The  viscous  aspects  of  transonic  flow  would  require 
Reynolds  numbers  to  be  specified.  Wing-body  interference  studies  require  a fuselage 
to  be  defined;  at  least  one  additional  external-store/pylon  configuration  should  also 
be  considered.  Control  surface  studies  would  require  the  planform  description  in 
addition  to  the  airfoil  thickness  and  Reynolds  numbers;  a definition  of  gaps  and  seals 
might  also  be  considered.  Rotary  effects  on  T-tails  can  be  calculated  for  the  existing 
standard  AGARD  configuration,  and  calculations  of  vortex  shedding  characteristics  can 
also  be  made  on  the  existing  standard  configurations.  Finally,  studies  of  rotor  blade 
loading  would  require  a standard  AGARD  rotor  blade  along  with  its  tip  Mach  number, 
advance  ratio,  and  related  parameters. 

It  is  hoped  that  some  of  the  above  topics  can  become  the  basis  of  a third  AGARD 
comparative  study  of  computed  results.  The  reference  list  presented  here  is  only  meant 
to  be  representative  of  the  state-of-the-art;  the  author  apologizes  to  the  Panel  for 
its  American  flavor. 
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Table  1 


Review  of  Calculation  Methods  for  Unsteady  Transonic  Flow 


Author 

Year 

Ref. 

Remarks 

Two-dimensional 

methods 

Stahara  and  Spreiter 

Isogai 

Dowell 

Revell 

Ehlers 

Beam  and  Warming 
Ballhaus  and  Lomax 
Chan  and  Brashears 
Traci,  Farr,  Albano  & Cheng 
Magnus  and  Yoshihara 
Isom  and  Caradonna 

1973 

1974 

1975 

1973 

1973 

1974 
1974 
1974 

1974 

1975 
1975 

15 

16,17 

18 

19 

20 
21 
22 

23 

24 

25,26 

27,28 

(M  as  1.0 
<no  shock 
(waves 

layered-medium 
finite  difference 
finite  difference 
finite  difference 
finite  element 
finite  difference 
finite  difference 
finite  difference 

Ruo  and  Theisen 

1973 

29 

(m  as  1.0 

Isogai 

1974 

17 

) no  shock  waves 

- 

— 

Three-dimensional 

Cunningham 

1973 

30-32 

mixed  subsonic-super- 
sonic method 

methods 

Garner 

1975 

33 

semi-empirical 

Isom  and  Caradonna 

1975 

27,28 

finite  difference 

Weatherhill,  Ehlers  and 
Sebastian 

1975 

34 

finite  difference 
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Aeroelasticity,  Vol.VI.  The  latter  effort  and  report  has  also  proved  to  be  highly  successful. 

^New  developments  are  rapidly  emerging  in  unsteady  aerodynamics.  The  aeroelastician  will 
continue  heavy  reliance  on  prediction  of  airloads  from  theoretical  methods.  A terse  description 
of  the  new  state  of  the  art  was  required  and  has  beerv  very  competently-provided  hy- 
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